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We derive an approximate expression for the Gilbert damping coefficient ac of itinerant electron 
ferromagnets which is based on their description in terms of spin-density-functional-theory (SDFT) 
and Kohn-Sham quasiparticle orbitals. We argue for an expression in which the coupling of mag- 
netization fluctuations to particle-hole transitions is weighted by the spin-dependent part of the 
theory's exchange-correlation potential, a quantity which has large spatial variations on an atomic 
length scale. Our SDFT result for ac is closely related to the previously proposed spin-torque 
correlation-function expression. 
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I. INTRODUCTION 

The Gilbert parameter ac characterizes the damping 
of collective magnetization dynamics^. The key role of 
ac in current-driven^ and precessional^ magnetization 
reversal has renewed interest in the microscopic physics 
of this important material parameter. It is generally 
accepted that in metals the damping of magnetization 
dynamics is dominated^ by particle-hole pair excitation 
processes. The main ideas which arise in the theory of 
Gilbert damping have been in place for some timo^iSi. It 
has however been difficult to apply them to real materi- 
als with the precision required for confident predictions 
which would allow theory to play a larger role in design- 
ing materials with desired damping strengths. Progress 
has recently been achieved in various directions, both 
through studies^ of simple models for which the damp- 
ing can be evaluated exactly and through analyses^ of 
transition metal ferromagnets that arc based on realis- 
tic electronic structure calculations. Evaluation of the 
torque correlation formula^ for ac used in the later calcu- 
lations requires knowledge only of a ferromagnet's mean- 
field electronic structure and of its Bloch state lifetime, 
which makes this approach practical. 

Realistic ah initio theories normally employ spin- 
density-functional theory^ which has a mean-field theory 
structure. In this article we use time-dependent spin- 
density functional theory to derive an explicit expression 
for the Gilbert damping coefficient in terms of Kohn- 
Sham theory eigenvalues and eigenvectors. Our final 
result is essentially equivalent to the torque-correlation 
formulai^ for ac, but has the advantages that its deriva- 
tion is fully consistent with density functional theory, 
that it allows for a consistent microscopic treatments of 
both dissipative and reactive coefficients in the Landau- 
Liftshitz Gilbert (LLG) equations, and that it helps 
establish relationships between different theoretical ap- 
proaches to the microscopic theory of magnetization 
damping. 

Our paper is organized as follows. In Section II 
we relate the Gilbert damping parameter uq of a fer- 
romagnet to the low-frequency limit of its transverse 



spin response function. Since ferromagnetism is due 
to electron-electron interactions, theories of magnetism 
are always many-electron theories, and it is necessary to 
evaluate the many-electron response function. In time- 
dependent spin-density functional theory the transverse 
response function is calculated using a time-dependent 
self- consistent-field calculation in which quasiparticles 
respond both to external potentials and to changes in the 
interaction-induced effective potential. In Section III we 
use perturbation theory and time-dependent mean-field 
theory to express the coefficients which appear in the 
LLG equations in terms of the Kohn-Sham eigenstates 
and eigenvalues of the ferromagnet's ground state. These 
formal expressions are valid for arbitrary spin-orbit cou- 
pling, arbitrary atomic length scale spin-dependent and 
scalar potentials, and arbitrary disorder. By treating dis- 
order approximately, in Section IV wc derive and com- 
pare two commonly used formulas for Gilbert damping. 
Finally, in Section V we summarize our results. 



II. MANY-BODY TRANSVERSE RESPONSE 
FUNCTION AND THE GILBERT DAMPING 
PARAMETER 

The Gilbert damping parameter ac appears in the 
Landau-Liftshitz- Gilbert expression for the collective 
magnetization dynamics of a ferromagnet: 



(1) 



In Eq.( [T|) Heff is an effective magnetic field which 
we comment on further below and Q = (QxT^y^^z) is 
the direction of the magnetization. This equation de- 
scribes the slow dynamics of smooth magnetization tex- 
tures and is formally the first term in an expansion in 
time-derivatives. 

The damping parameter ac can be measured by per- 
forming ferromagnetic resonance (FMR) experiments in 
which the magnetization direction is driven weakly away 
from an easy direction (which we take to be the z- 
dircction.). To relate this phenomenological expression 



2 



formally to microscopic theory we consider a system in 
which external magnetic fields couple onlyii to the elec- 
tronic spin degree of freedom and associate the magneti- 
zation direction n with the direction of the total electron 
spin. For small deviations from the easy direction, Eq.([T]) 
reads 
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The gyromagnetic ratio has been absorbed into the units 



of the field 'Heff so that this quantity has energy units 
and we set h = 1 throughout. The corresponding formal 
linear response theory expression is an expansion of the 
long wavelength transverse total spin response function 
to first order— in frequency uj: 
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where a, /? e {x, y}, lu = idt is the frequency, 5*0 is the to- 
tal spin of the fcrromagnet, Tiext is the external magnetic 
field and x is the transverse spin-spin response function: 
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Here \^n) is an exact eigenstate of the many-body Hamiltonian and uJn,o is the excitation energy for state n. We use 
this formal expression below to make some general comments about the microscopic theory of ac- In Eq.([3l) Xa0 
the static (w = 0) limit of the response function, and ^ is the first derivative with respect to uj evaluated at a; = 0. 
Notice that we have chosen the normalization in which x is the total spin response to a transverse field; x is therefore 
extensive. I 



The key step in obtaining the Landau-Liftshitz-Gilbert 
form for the magnetization dynamics is to recognize that 
in the static limit the transverse magnetization responds 
to an external magnetic field by adjusting orientation to 
minimize the total energy including the internal energy 
Eint and the energy due to coupling with the external 
magnetic field, 



It follows that 
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We obtain a formal equation for H^ff corresponding to 
Eq.([2) by multiplying Eq.([31) on the left by [Xa fj]~^ ^^'^ 
recognizing 
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as the internal energy contribution to the effective mag- 
netic field "Heff = 'Hint + Hext- With this identification 
Eq.([3]) can be written in the form 



H, 



eff. 
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where 



Ca,p - -SoHx'T' X' (x^*)-']a,/3 = lSod^Xa!l3- (9) 



According to the Landau-Liftshitz Gilbert equation then 
r 

*^x,y 



-Cy^x = 1 and 
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Explicit evaluation of the off-diagonal components of C 
will in general yield very small deviation from the unit 
result assumed by the Landau-Liftshitz-Gilbert formula. 
The deviation refiects mainly the fact that the magneti- 
zation magnitude varies slightly with orientation. We do 
not comment further on this point because it is of little 
consequence. Similarly Cx,x is not in general identical 
to Cy^y, although the difference is rarely large or impor- 
tant. Eq.f fTU]) is the starting point we use later to derive 
approximate expressions for uq. 

In Eq.([9]) Xa,i3{^) is the correlation fimction for an 
interacting electron system with arbitrary disorder and 
arbitrary spin-orbit coupling. In the absence of spin- 
orbit coupling, but still with arbitrary spin-independent 
periodic and disorder potentials, the ground state of a 
ferromagnet is coupled by the total spin-operator only to 
states in the same total spin multiplet. In this case it 
follows from Eq. ( |4]) that 
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where Hq is a static external field, which is necessary 
in the absence of spin-orbit coupling to pin the magne- 
tization to the z direction and splits the fcrromagnet 's 
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ground state many-body spin multiplet. Similarly 



exchange-correlation potential and 
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1, yielding 
Spin-orbit 



coupling is required for magnetization damping^. 



III. SDF-STONER THEORY EXPRESSION FOR 
GILBERT DAMPING 

Approximate formulas for ac in metals are inevitably 
based on on a self-consistent mean-field theory (Stoner) 
description of the magnetic state. Our goal is to derive an 
approximate expression for uq when the adiabatic local 
spin-density approximation^, is used for the exchange cor- 
relation potential in spin-density-functional theory. The 
effective Hamiltonian which describes the Kohn-Sham 
quasiparticle dynamics therefore has the form 



Uks = np- A(n(f), \s{r)\) n{f) ■ s, 



(13) 



where Tip is the Kohn-Sham Hamiltonian of a paramag- 
netic state in which |s(r)|(the local spin density) is set to 
zero, s is the spin-operator, and 



A(n, s) 



d [nexc{n,s)] 
ds 
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is the magnitude of the spin-dependent part of the 
exchange-correlation potential. In Eq.f e2;c(n, s) is 
the exchange-correlation energy per particle in a uni- 
form electron gas with density n and spin-density s. 
We assume that the ferromagnet is described using 
some semi-relativistic approximation to the Dirac equa- 
tion like those commonly usedi^ to describe magnetic 
anisotropy or XMCD, even though these approximations 
are not strictly consistent with spin-density-functional 
theory. Within this framework electrons carry only a 
two-component spin-1/2 degree of freedom and spin-orbit 
coupling terms are included in Tip. Since nea;c(n, s) ~ 
[(n/2 + s)^/^ + {n/2 - sY'^], An{n,s) ~ n^/^ is larger 
closer to atomic centers and far from spatially uniform on 
atomic length scales. This property figures prominently 
in the considerations explained below. 

In SDFT the transverse spin-response function is ex- 
pressed in terms of Kohn-Sham quasiparticle response to 
both external and induced magnetic fields: 



dr' 



so(r)f]„(f) = J ^QP{r,r') [HextMr')+MrO^0ir')]. 

(15) 

In Eg.l' lisp V is the system volume, so(^) is the magni- 
tude of the ground state spin density, Ao(r) is the mag- 
nitude of the spin-dependent part of the ground state 
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where fi is the ground state Kohn-Sham occupation fac- 
tor for eigenspinor |i) and ojij = — ej is a Kohn- 
Sham eigenvalue difference, ^'^^(r, r') has been normal- 
ized so that it returns the spin-density rather than total 
spin. Like the Landau-Liftshitz-Gilbert equation itself, 
Eq.( (TS]) assumes that only the direction of the mag- 
netization, and not the magnitudes of the charge and 
spin-densities, varies in the course of smooth collective 
magnetization dynamicsi^. This property should hold 
accurately as long as magnetic anisotropics and exter- 
nal fields are weak compared to Aq. We are able to use 
this property to avoid solving the position-space integral 
equation imphed by Eq.([T5]). Multiplying by Ao(r) on 
both sides and integrating over position we findi^ that 
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where we have taken advantage of the fact that in FMR 
experiments TLext,i3 and are uniform. Aq is a spin- 
density weighted average of Ao(7"), 
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is the response function of the transverse-part of the 
quasiparticle exchange-correlation effective field response 
function, not the transverse-part of the quasiparticle 
spin response function. In Eq.( [TO]) . {i\0{f)\j) = 
J dfO{f){i\r){f\j) denotes a single-particle matrix ele- 
ment. Solving Eq.f I17p for the many-particle transverse 
susceptibility (the ratio of So^a to Hext,f3) and inserting 
the result in Eq. ( [5]) yields 



(20) 



Our derivation of the LLG equation has the advantage 
that the equation's reactive and dissipative components 
are considered simultaneously. Comparing Eq. ( llSp and 
Eq.([7]) we find that the internal anisotropy field can also 
be expressed in terms of x*^^- 



E[x?S- 



^(c. = 0)- 



ftp. (21) 



Eq.( 1^01) and Eq.( [?T|) provide microscopic expressions 
for all ingredients that appear in the LLG equations 
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linearized for small transverse excursions. It is gener- 
ally assumed that the damping coefficient ac is inde- 
pendent of orientation; if so, the present derivation is 
sufficient. The anisotropy- field at large transverse ex- 
cursions normally requires additional information about 
magnetic anisotropy. We remark that if the Hamiltonian 
does not include a spin-dependent mean-field dipole in- 
teraction term, as is usually the case, the above quantity 
will return only the magnetocrystalline anisotropy field. 
Since the magnetostatic contribution to anisotropy is al- 
ways well described by mean-field-theory it can be added 
separately. 

We conclude this section by demonstrating that the 
Stoner theory equations proposed here recover the exact 
results mentioned at the end of the previous section for 
the limit in which spin-orbit coupling is neglected. We 
consider a SDF theory ferromagnet with arbitrary scalar 
and spin-dependent effective potentials. Since the spin- 
dependent part of the exchange correlation potential is 
then the only spin-dependent term in the Hamiltonian it 
follows that 

[H-KS, Sa] = f^a.p Ao(r)s;3 (22) 
and hence that 

{i\saAo{r)\j) = -~iea,pUJtj{i\si3\j). (23) 

Inserting Eq.( I23p in one of the matrix elements of 
Eq. ( [TO]) yields for the no-spin-orbit-scattering case 

xZp(''^0)^Sa.,p SoAo. (24) 



The internal magnetic field 'Hint,a is therefore identically 
zero in the absence of spin-orbit coupling and only exter- 
nal magnetic fields will yield a finite collective precession 
frequency. Inserting Eq. ( I23p in both matrix elements of 
Eq.([T9l) yields 

du,Mx^2 = ea,/35o. (25) 

Using both Eq. ( [24]) and Eq. ( [25l) to invert x'^^ we re- 
cover the results proved previously for the no-spin-orbit 
case using a many-body argument: Cx,y = '^^y.x ~ 1 
and Cx.x = ^y.y = 0. The Stoner-theory equations de- 
rived here allow spin-orbit interactions, and hence mag- 
netic anisotropy and Gilbert damping, to be calculated 
consistently from the same quasiparticle response func- 
tion x*^^- 



IV. DISCUSSION 

As long as magnetic anisotropy and external magnetic 
fields arc weak compared to the exchange-correlation 
splitting in the ferromagnet we can use Eq. ( I24|l to ap- 
proximate Xa^i^ ~ Using this approximation and 
assuming that damping is isotropic we obtain the follow- 
ing explicit expression for temperature T — > 0: 



ao = Cx,x = -SoAld^lmlx^^-^] = ^ V 6{ej ~ ep) <5(e, - ep) (j|s.Ao(r)|i) {i\sxAo{f)\j) 

= f E ^('1 - ^('^ - {mp.^vm (imp,sy]\j). (26) 

•JO 

The second form for aa is equivalent to the first and follows from the observation that for matrix elements between 
states that have the same energy 

(i|[WifS,Sa]|j) = -iea,/3 (i|Ao(r)s/3|i) + (i|[Hp,Sa]|j) = (for U^J = 0). (27) 

Eq. ( I26p is valid for any scalar and any spin-dependent potential. It is clear however that the numerical value of aa 
in a metal is very sensitive to the degree of disorder in its lattice. To see this we observe that for a perfect crystal 
the Kohn-Sham eigenstates are Bloch states. Since the operator Ao(r)sQ, has the periodicity of the crystal its matrix 
elements are non-zero only between states with the same Bloch wavevector label k. For the case of a perfect crystal 
then 

^-f) '^('^fcn ~ ''f) {kn'\sxAii{r)\kn) (fcn|s^Ao(r)|fcn') 

ep) '5(e£„ -ef) {kn'\[np,Sy\\kn) {kn\[np, Sy\\kn') . (28) 

I 

where nn' are band labels and sg is the ground state Brillouin-zonc (BZ). 
spin per unit volume and the integral over k is over the 
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Clearly ac divergesi^ in a perfect crystal since 
(fcn|Sa;Ao(r)|fcn) is generically non-zero. A theory of 
ac rnust therefore always account for disorder in a crys- 
tal. The easiest way to account for disorder is to replace 
the ^(e^„ — e_F) spectral function of a Bloch state by a 
broadened spectral function evaluated at the Fermi en- 
ergy A^J^tp). If disorder is treated perturbatively this 
simple ansatz can be augmentedi^ by introducing impu- 
rity vertex corrections in Eq. ( [^5)1 . Provided that the 
quasiparticle lifetime is computed via Fermi's golden rule, 



these vertex corrections restore Ward identities and yield 
an exact treatment of disorder in the limit of dilute im- 
purities. Nevertheless, this approach is rarely practical 
outside the realm of toy models, because the sources of 
disorder are rarely known with sufficient precision. 
Although appealing in its simplicity, the <5(ej:„ — e^?) — > 
Aji^^ep) substitution is prone to ambiguity because it 
gives rise to qualitatively different outcomes depending 
on whether it is applied to the first or second line of Eq. 
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So Jbz (271")^ 
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E ^£,n {£F)Aj:^^,ieF){kn'\['Hp,Sy]\kn){kn\[np,.Sy]\ kn') , 
^ (ej. (ej. ) (fcn' I Ao (r) I fcn) (fcn I Ao (f ) I fcn'} . 
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ajj is the torque-correlation (TC) formula used in 

realistic electronic structure calculations^ and cSq^^ is 
the spin-flip (SF) formula used in certain toy model 
calculationsi^. The discrepancy between TC and SF ex- 
pressions stems from inter-band (n ^ n') contributions 
to damping, which may now connect states with dif- 
ferent band energies due to the disorder broadening of 
the spectral functions. Therefore, {kn\\HKSiSa\\kn') no 
longer vanishes for n ^ n' and Eq. ( [27|) indicates that 
^.C^C") ^ a^Q^^ only if the Gilbert damping is dominated 
by intra-band contributions and/or if the energy differ- 
ence between the states connected by inter-band transi- 
tions is small compared to Aq. When Oq'^^ ^ ckcf^\ 
it is a priori unclear which approach is the most accu- 
rate. One obvious flaw of the SF formula is that it pro- 
duces a spurious damping in absence of spin-orbit interac- 
tions; this unphysical contribution originates from inter- 
band transitions and may be cancelled out by adding 
the leading order impurity vertex corrections^. In con- 
trast, \Hp,Sy\ = in absence of spin-orbit interaction 
and hence the TC formula vanishes identically, even with- 
out vertex corrections. From this analysis, TC appears 
to have a pragmatic edge over SF in materials with weak 
spin-orbit interaction. However, insofar as it allows inter- 
band transitions that connect states with w^.j > Aq, 
TC is not quantitatively reliable. Furthermore, it can 



be shownii that when the intrinsic spin-orbit coupling 
is significant (e.g. in ferromagnetic semiconductors), the 
advantage of TC over SF (or vice versa) is marginal, and 
impurity vertex corrections play a significant role. 

V. CONCLUSIONS 

Using spin-density functional theory we have derived 
a Stoner model expression for the Gilbert damping co- 
efficient in itinerant ferromagnets. This expression ac- 
counts for atomic scale variations of the exchange self 
energy, as well as for arbitrary disorder and spin-orbit 
interaction. By treating disorder approximately, we have 
derived the spin-flip and torque-correlation formulas pre- 
viously used in toy-model and ab-initio calculations, re- 
spectively. We have traced the discrepancy between these 
equations to the treatment of inter-band transitions that 
connect states which are not close in energy. A better 
treatment of disorder, which requires the inclusion of im- 
purity vertex corrections, will be the ultimate judge on 
the relative reliability of either approach. When damping 
is dominated by intra-band transitions, a circumstance 
which we believe is common, the two formulas are identi- 
cal and both are likely to provide reliable estimates. This 
work was suported by the National Science Foundation 
under grant DMR-0547875. 



For a historical perspective see T.L. Gilbert, IEEE Trans. ^ J.A.C. Bland and B. Heinrich (Eds.), Ultrathin Mag- 

Magn. 40, 3443 (2004). netic Structures III: Fundamentals of Nanomagnetism 

^ For an introductory review see D.C. Ralph and M.D. Stiles, (Springer- Verlag, New York, 2005). 

J. Magn. Mag. Mater. 320, 1190 (2008). * V. Korenman and R. E. Prange, Phys. Rev. B 6, 2769 



6 



(1972). 

^ V. Kambersky, Czech J. Phys. B 26, 1366 (1976). 

^ Y. Tserkovnyak, G.A. Fiete, and B.I. Halperin, Appl. 
Phys. Lett. 84, 5234 (2004); E.M. Hankiewicz, G. Vig- 
nale and Y. Tserkovnyak, Phys. Rev. B 75, 174434 (2007); 
Y. Tserkovnyak et al, Phys. Rev. B 74, 144405 (2006) ; 
H.J. Skadsem, Y. Tserkovnyak, A. Brataas, G.E.W. Bauer, 
Phys. Rev. B 75, 094416 (2007); H. Kohno, G. Tatara 
and J. Shibata, J. Phys. See. Japan 75, 113706 (2006); 
R.A. Duine et al, Phys. Rev. B 75, 214420 (2007). Y. 
Tserkovnyak, A. Brataas, and G.E.W. Bauer, J. Magn. 
Mag. Mater. 320, 1282 (2008). 

K. Gilmore, Y.U. Idzerda and M.D. Stiles, Phys. Rev. Lett. 
99, 27204 (2007); V. Kambersky, Phys. Rev. B 76, 134416 
(2007). 

® For zero spin-orbit coupling aa vanishes even in presence 
of magnetic impurities, provided that their spins foUow the 
dynamics of the magnetization adiabatically. 

^ O. Gunnarsson, J. Phys. F 6, 587 (1976). 
1° Z. Qian, G. Vignale, Phys. Rev. Lett. 88, 056404 (2002). 
In doing so we dodge the subtle difficulties which compli- 
cate theories of orbital magnetism in bulk metals. See for 
example J. Shi, G. Vignale, D. Xiao, and Q. Niu, Phys. 
Rev. Lett. 99, 197202 (2007); I. Souza and D. Vanderbilt, 
Phys. Rev. B 77, 054438 (2008) and work cited therein. 
This simplification should have little influence on the the- 
ory of damping because the orbital contribution to the 
magnetization is relatively small in systems of interest and 



because it in any event tends to be coUinear with the spin 
magnetization. 

For most materials the FMR frequency is by far the small- 
est energy scale in the problem. Expansion to linear order 
is almost always appropriate. 

See for example A.C. Jenkins and W.M. Temmerman, 
Phys. Rev. B 60, 10233 (1999) and work cited therein. 
This approximation does not preclude strong spatial varia- 
tions of |so(f01 ^^'^ l^o(r)l atomic lenghtscales; rather it 
is assumed that such spatial profiles will remain unchanged 
in the course of the magnetization dynamics. 
For notational simplicity we assume that all magnetic 
atoms are identical. Generalizations to magnetic com- 
pounds are straight forward. 

Eq. ([26} is valid provided that cut << 1. While this con- 
dition is normally satisfied in cases of practical interest, it 
invariably breaks down as r — > oo. Hence the divergence 
of Eq. ( I26|l in perfect crystals is spurious. 
I. Garate and A.H. MacDonald (in preparation). 
J. Sinova et al, Phys. Rev. B 69, 85209 (2004). In order to 
get the equivalence, trade hz by Ao and use Ao = JpdSo, 
where Jpd is the p-d exchange coupling between GaAs va- 
lence band holes and Mn d-orbitals. In addition, note that 
our spectral function differs from theirs by a factor 2-k. 
H. Kohno, G. Tatara and J. Shibata, J. Phys. Soc. Japan 
75, 113706 (2006). 



